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Abstract. We have constrained the extended (delayed and accelerated) models of 
hydrogen recombination, by investigating associated changes of the position and the 
width of the last scattering surface. Using the recent CMB and SDSS data, we find that 
the recent data constraints favor the accelerated recombination model, though the other 
models (standard, delayed recombination) are not ruled out at 1-cr confidence level. If 
the accelerated recombination had actually occurred in our early Universe, baryonic 
clustering on small-scales is likely to be the cause of it. By comparing the ionization 
history of baryonic cloud models with that of the best-fit accelerated recombination 
model, we find that some portion of our early Universe has baryonic undcrdensity. We 
have made the forecast on the PLANCK data constraint, which shows that we will be 
able to rule out the standard or delayed recombination models, if the recombination in 
our early Universe had proceeded with ~ —0.01 or lower, and residual foregrounds 
and systematic effects are negligible. 
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1. Introduction 

The recombination process of cosmic plasma, which have occurred around the redshift 
Zrec — 1100, decouples photons from baryons [H l2]- In the presence of Ly-a 
photon sources, the recombination process might proceed with delay [21 HI [5], or with 
acceleration in the presence of baryonic clustering on small-scales [HI [71 Ej • The Cosmic 
Microwave Background (CMB) anisotropy, which is sensitive to the ionization history 
of our Universe, is affected by delay or acceleration of the recombination process 
[31 m El [6l [7j. The major effect on the CMB anisotropy is expected in the CMB 
polarization and small angular scale temperature anisotropy [H El El [8] . The five year 
data of the Wilkinson Microwave Anisotropy Probe (WMAP) [9l [TOl [11] is released 
and the recent ground-based CMB observations such as the ACBAR [121 113] and 
QUaD [m [151 [IB] provide information complementary to the WMAP data. In near 
future, PLANCK surveyor [17] is going to measure CMB temperature and polarization 
anisotropy with great accuracy over wide range of angular scales. In this paper, we 
investigate the recombination models and baryonic cloud models, using the recent CMB 
and SDSS data. 

The outline of this paper is as follows. In Section [21 we investigate the 
extended recombination models and shows that the data constraints favor accelerated 
recombination models. In Section [3l we discuss and constrain baryonic cloud models. 
In Section [H we make the forecast on the PLANCK data constraint. In Section [5l we 
summarize our investigation with conclusion. 

2. Distortion on the standard ionization history 

The presence of extra resonance photon sources may delay the recombination of cosmic 
plasma [31 [H [5l [6], while the presence of baryonic clustering on small scales may 
accelerate the recombination [SI [71 [8] . The simplest model for the production of extra 
resonance photons is given by [31 [6] : 

^ = e^{z)H{z)n, (1) 

where n is the number density of atoms, H{z) is the Hubble expansion rate at a redshift 
z, and ea{z) is a parameter dependent on the production mechanism. Since thickness of 
the last scattering surface is very small in comparison to the horizon of the last scattering 
surface Li^, the dependence of ea{z) on z can be parametrized as ea{zrec) + o{A/Lis). 
Hence, we use the approximation ea{z) ^ ea{zrec) through our investigation. Though 
Eq. [T] was originally proposed for the delayed recombination, we use Eq. [T] to model 
the accelerated recombination by assigning negative values to [6l [71 [8] . However, it 
should be noted that the physical basis of the accelerated recombination (i.e. baryonic 
clustering on small scales) is different from that of the delayed recombination. We 
also would like to point out that the distortion on the CMB black body spectrum 
by extended recombination process (0 < Icdl < 1) is negligible in comparison with the 
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distortion by the re-ionization (see [18j for details), and is well within the COBE FIRAS 
data constraint [19]. 
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Figure 1. Ionization history for ea = 0.3, 0.1, 0, —0.03, —0.07 (from the highest curve 
to the lowest) 
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Figure 2. Temperature anisotropy power spectrum for = —0.07, —0.03, 0, 0.1, 0.3 
in the descending order of Doppler peak heights. The ACBAR 2008 and the binned 
WMAP 5 year data are shown with error bars. 



By making a small modification to the RECFAST code [20l [HI [22] , we have computed 
the ionization fraction Xe for various and plotted them in Fig. [H We may see that 
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the ionization fraction Xg of the recombination models Cq, < drops much faster than 
the standard model ea = 0, while Xe of the recombination models > drops much 
slower. 

150^ • 
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Figure 3. Temperature-E mode correlation for = —0.07, —0.03, 0, 0.1, 0.3 in the 
ascending order of trough heights. The QUaD data are shown with error bars. 
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Figure 4. E mode power spectrum for = —0.07, —0.03, 0, 0.1, 0.3 in the 
descending order of Doppler peak heights. The QUaD data are shown with error 
bars. 



In Fig. [21 [3] and HJ we show the temperature power spectra, E mode power spectra 
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and TE correlation for various e^- It is worth noting that the location and heights of 
the Doppler peaks are affected by the delay and acceleration of the recombination. 
As shown in Fig. [21 the accelerated recombination models {ea < —0.07) are not 
in good agreement with the WMAP data constraint (e.g. the first Doppler peak), 
and the delayed recombination models (e^ > 0.3) are ruled out PJ. Hence we have 
investigated the recombination models (—0.07 < < 0.3). For data constraints, 
we have used the Sloan Digital Sky Survey (SDSS) data [23l [2ll [25] and the recent 
CMB observations (the WMAP 5 year data P HD], the ACBAR 2008 [El [I3] and the 
QUaD [T5l [16] ). Through small modifications to the CosmoMC package [26], we have 
included the parameter ea of the prior distribution — 0.07<ea<0.3in the cosmo logical 
parameter estimation and explored the multi-dimensional parameter space {^Ibh"^, VLji^, 
T, Us, log[10^°y4s], Hq, ea) by fitting the matter power spectra to the SDSS data, and 
the CMB anisotropy power spectra Cf^, Cf^ and Cf^ to the recent CMB observations 
(WMAP5YR + ACBAR + QUaD). 




Figure 5. CMB + SDSS : the marginalized likelihood (the solid curve), the mean 
likelihood (the dashed curve). The normalization is chosen such that the peak value is 
equal to unity. 



We have run the modified CosmoMC on a MPI cluster with 6 chains. For the 
convergence criterion, we have adopted the Gelman and Rubin's "variance of chain 
means" and set the R-1 statistic to 0.03 for stopping criterion [271 [2B]- The convergence, 
which is measured by the R-1 statistic, is 0.0262 and 15107 chains steps are used. 

In Fig. [3 we show the marginalized likelihood and mean likelihood of e^, given 
those observations. The solid lines and dotted lines correspond to the marginalized 
likelihood and the mean likelihood respectively (for distinction between marginalized 
likelihood and mean likelihood, refer to [26]). Though we are aware that CosmoMC 
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Figure 6. The marginalized likelihood in the plane of Ca versus 6 basic parameters, 
using the CMB + SDSS constraint. Two contour lines correspond to la and 2a levels. 



does not provide a precise best-fit value, we quote the best-fit value from CosmoMC 
as often done in literature. Using the recent CMB + SDSS observation constraints, 
we find ea = -0.003421°;°^^^ at la confidence level and e„ = -0.00342+°;°^^^ at 2a 
confidence level. As also shown in Fig. [5l the current data constraints favor accelerated 
recombination models, though eo, > is still in la confidence interval. In Fig. [6], we 
have plotted the marginalized distribution and mean likelihoods in the plane of versus 
{Qbh"^, ^ch'^, T, Ug, \og[10^^ As], Hq}. We summarize la constraints on cosmological 
parameters in Table [H given CMB + SDSS data. In comparison to the cosmological 
parameters estimated with the standard recombination model [Hj, optical depth r is 
affected most, while other cosmological parameters are affected negligibly. 



Table 1. cosmological parameters of ACDM + the extended recombination model 
constrained by WMAP5YR + ACBAR + QUaD + SDSS data 

parameter la constraint 

Vtbh^ 0.0229l°-°°}7 



T 0.0930.0.0448 



+0.0384 
0.0448 
0.0316 



n, 0.9625_o.o3i7 
log[W^As] at ko = 0.005/Mpc 3.0547+°;°^^^ 

Ho 71.73031^;^^^^ km/s/Mpc 
. _n nnQ4+0-0i85 

ta U.UUO^_o.0237 
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3. Small-scale baryonic cloud models and accelerated recombination 

As discussed in |7], tlie presence of the baryonic clustering in the range of very small 
mass scales M ~ 10 — IO^Mq can cause accelerated recombination. The simplest model 
for baryonic clustering is given by a baryonic cloud model [7] as follows: 

Pb = Pb,inf + Pb,out{l - f) (2) 

, where is the mean density of baryonic matter, and Ph,out and pfe^j„ are the baryonic 
density of baryonic clouds and intercloud regions respectively, and / is the total 
volume fraction of intercloud regions. Denoting the baryoninc density contrast between 
intercloud region and clouds by ^ = Pb,in/ Pb,out, it may be easily shown that 

- ^Pb _ Pb (o\ 

Pb,in ^^f^^_^y Pb,out ^_^f^^_^y 

On the other hand, in the presence of baryonic clouds, there arises diffusion from clouds 
into intercloud regions, washing out density contrast. The characteristic scale of such 
diffusion is close to the Jean's length, Rj ~ Cgrjr, where is the baryonic speed of sound 
and rjr = J cdt/a{t) is the conformal time corresponding to the recombination time. 
Baryonic clouds of mass scales M ~ 10 — IO^Mq has length scale Ra < Rj < R < Ry, 
where Ra and R^ are the mean free path of resonance photons and CMB photons 
respectively. Hence we may neglect baryonic diffusion. Free electron density of clouds 
and intercloud regions, which we denote by Ue^out and ne,m respectively, are given by 

- n,,„, n.,., - ^1 - n,,^„ (4J 

where Xe^out and Xe^m are the ionization fraction of clouds and intercloud regions, and 
Yp and rib are the mass fraction of ^He and baryon number density respectively. The 
ionization fraction Xe is given by ne/{ne + ^p), where Up is the total number density of 
free protons and protons trapped in nucleus. Since scales of baryonic clouds are smaller 
than the mean free path of CMB photons, the effective ionization fraction, which CMB 
anisotropy is sensitive to, is given by the mean ionization fraction: 



where the mean value of free electron density is given by [7]: 

(rie) = rie^inf + ^^e,o«^(l - /)• 

Using Eq. [2], H] and Eq. [5l we may show that the effective ionization fraction has the 
following relation to the ionization fraction of clouds and intercloud regions: 

{Xe} Xe^inGin ~\~ Xf^^outG outi (6) 



where 



G.n = , , i{ ( \ !^^T(f ) , (7) 



l + /(e-l) \l-{Yp)/2 
n — ^ ^ f ( 1 ~ Yp^out/'^ \ 

l + /(e-l) Vl-(^p)/2 ^ ^ 
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We would like to remind that the effective ionization fraction (xg) is not necessarily 
equal to the ionization fraction of homogeneous baryonic model (i.e. ^ = 1), because 
the recombination process has non-linear dependence on baryonic density. 




1+z 

Figure 7. Ionization history of the accelerated recombination model (eq = —0.0034) 
and the baryonic cloud model (/ — 0.019, ^ = 0.038): the solid curve shows the 
ionization history of the accelerated recombination model and the dashed curve shows 
ionization history of a baryonic cloud model. Two curves are visually identical. 

As presented in the previous section, the CMB data constraints with or without 
SDSS data favor the accelerated recombination models. Given certain values of / and 
^, we can compute ionization history (i.e. Xe{z)) of cloud and intercloud region, using 
the RECFAST with the respective values of Yp and fib for cloud and intercloud region. 
Exploring two dimensional parameter space (/, ^), we have fitted the ionization history 
of baryonic cloud models to the best-fit accelerated recombination model. In Fig. [7l we 
show the ionization fraction Xe as a function of redshift z. The solid curve shows the 
ionization fraction Xg of the accelerated recombination model (e^ = —0.0034), while the 
dashed curve shows the mean ionization fraction of the fitted baryonic cloud model 
(/ = 0.019, ^ = 0.038). As shown in Fig. [3, two curves in Fig. [7] are visually 
indistinguishable and the fitting errors are within the accuracy of RECFAST [29j. 



Table 2. Best- fit baryonic cloud model 



Constraints 


ea / 




Pb,in/ Pb 


Pb,out/ Pb 


CMB + SDSS 


-0.0034 0.019 


0.038 


0.04 


1.02 



Provided that = —0.0034, we find that baryonic clouds with baryonic density 
1.02 Pb occupy ~ 98% of the total volume in our early Universe, while intercloud regions 
with baryonic density 0.04 pb occupy ~ 2%. In Table. [21 we summarize the best-fit 
values of /, ^, Pb,in/Pb and pb,out/Pb- 
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4. Forecast on the PLANCK data constraint 

In this section, we forecast the PLANCK data constraint on ea- We have assumed 
FWHM=10', AT/T = 2.5 (Stokes I) and AT/T = 4 (Stokes Q&U) [HI |30]. The pixel 
size is assumed to be equivalent to HEALPix [3T1 [32] pixellization of Nside=1024 (Res 
10). Assuming isotropic noise, we compute noise power spectrum as follows: 

A^z = iVoAfie''"' (9) 

, where a = FWHM//(8 ln2), A^'o is noise variance per pixel and Afi is the solid angle 
of a single pixel [Ij. 

We have obtained the forecast by feeding the PLANCK mock data and the SDSS 
data to the CosmoMC. The Planck mock data are generated by drawing spherical harmonic 
coefficients (2 < / < 1500) of signal and noise respectively from Gaussian distributions, 
whose variance are equal to the CMB power spectra and Eq. M respectively. We have 
obtained the CMB power spectra by using CAMB [33] with the modified RECFAST and 
ea = —0.01. In the Fig. [HI we show the marginalized likelihood (the solid curve) and 
the mean likelihoods (the dashed curve) of e^, given the Planck mock data. 




Figure 8. Planck mock data (eq = —0.01) : the marginalized likelihood (the solid 
curve), the mean likelihood (the dashed curve). The normalization is chosen such that 
the peak value is equal to unity. 

Given the mock data, the CosmoMC finds —0.015 < < —3.5 x 10"^ and 
—0.017 < ta < 0.0023 at la and 2a level respectively, which is similar to the forecast 
by [5j. While the CosmoMC forecast in Fig. [8] is made with the mock data of = —0.01, 
we need to make forecast for of various values. However, it is not practically feasible 
to run CosmoMC for various because of time usually required for running CosmoMC. 
Hence, we have estimated la interval for various via the Fisher matrix [ll[3ll, which 
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is given by [U El] : 



{- 



dXidXj 



N 



(10) 
(11) 



where Aj denotes the parameters to be estimated. Evaluated at the maximum of the 
hkehhood, the square root of diagonal element of the inverse Fisher matrix yields the 
marginalized 1 — a error on the parameter estimation [U [M]. The likelihood function 
for temperature and polarization anisotropy may be written as follows [35] : 



C 



{2n 



^ rexp[-iA(S + N)-^A], 

2 



where N is the number of data, A is a data vector, S is a signal covariance matrix 
and N is a noise covariance matrix. We expand anisotropy map in spherical harmonics 
so that our data vector A consists of spherical harmonic coefficients af^ and af^. In 
general, the computation of Eq. [TT] for the high resolution whole-sky observation such 
as the PLANCK observation is not possible on the modern computer at this moment. 
To facilitate the computation, we have made a few approximations. We assume very 
effective foreground cleaning with no need for sky masking and uniform instrument noise 
so that the signal and noise covariance matrices are block-diagonal and do no depend on 
spherical harmonic index m. In such configuration, signal and noise covariance matrices 
are given by: 



S = diag 



(jTT 




and 



N = diag 




v 




Taking into account the repeating pattern of diagonal blocks, we may show 
1 9S 



where 



dXi 



(S + N) 



-ii 



lTr[L#(S + N)-^|^fS 



dX, 



N) 



-ii 



(12) 



S = diag 









21 + 


.)) 


(jTT 
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N 



The right hand side of Eq. [12] can be easily computed in reasonable amount of time even 
for Post-PLANCK whole-sky observations. We chose six basic parameters (fift/i^, Qch"^, 
h, As, rig) plus for estimation parameters and set the values of basic six parameters 
to the WMAP best-fit values [H], [35l [36] . Due to non-linear dependence of CMB power 
spectra on parameter Aj, we resorted to numerical differentiation to obtain dS/dXi. The 
numerical derivatives are obtained by computing the following: 



q(\ + |AA,)-q(\-iAA,) 
AA,, 



(13) 



where we have set AAj/Aj = 10 ^ for six parameters. For the given set of parameter Aj, 
we have computed Cf^, Cf^, Cf ^, using CAMB ^ with the modified RECFAST. While 
we have made the forecast for various ea in the range of —0.07 < < —0.3, we find 
that Eq. [13] has numerical instability for Ae^ of very small values. Hence, we fixed Ae^ 
to be —0.001 instead of setting Aea/e^ = 10^'^. 
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Figure 9. forecast on estimation from the PLANCK data 



The left plot in Fig. [9] shows the estimation error for —0.07 < < 0.3, while the 
right one shows the estimation error in the vicinity of = 0. Two solid lines denote the 
boundary of Icr confidence interval and the dashed line shows the central values of la 
interval, which are set to the true values of ea- We find that la error tends to increase 
with decreasing and approaches ~ 0.005. We also find that the la error forecast 
by Fisher matrix method is smaller than the la forecast by CosmoMC. The discrepancy 
between CosmoMC forecast and Fisher matrix forecast is attributed to the deviation of 
parameter likelihood from Gaussian distribution. While Fisher matrix forecast is based 
on the marginalized likelihood with Gaussian approximation, CosmoMC forecast estimates 
la and 2a interval from mean likelihood 0. However, we can learn from the Fisher matrix 
forecast that la error for —0.7 < Cq, < 0.3 are at the same order of magnitude, which 



I For Gaussian distribution, marginalized likelihood and mean likelihood are identical. 
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provides complementary information to the CosmoMC forecast. For CosmoMC forecast as 
well as Fisher matrix forecast, we have neglected residual foregrounds and systematic 
effects including anisotropy of instrument noise and beam asymmetry, which will be 
additional effective sources of noise. Hence, our forecast should be regarded as more of 
a lower limit on the estimation variance. 

5. Discussion 

We have investigated the extended recombination models, using the recent CMB and 
SDSS data. We find that the data constraints favor accelerated recombination models, 
though other recombination models (standard, delayed recombination) are not ruled out 
at la confidence level. By comparing the ionization history of baryonic cloud models 
with the best-fit accelerated recombination model, we have constrained the baryonic 
cloud models, from which we find that our early Universe might have slight overdensity 
of baryonic matter ~ 1.02 pb for ~ 98% of total volume and underdensity of baryonic 
matter ~ 0.04 pb for the rest of space. The origin of primordial baryonic clouds, if exists, 
might be associated with inhomogeneous baryogenesis [37] in our very early Universe. 

While we have constrained baryonic cloud models indirectly by fitting the ionization 
history, more accurate constraint will be obtained only when CMB power spectra of the 
baryonic cloud models are fitted directly to the data. Once we get enough evidence 
of accelerated recombination from the upcoming PLANCK data, we plan to constrain 
baryonic clouds models directly. 
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